Conductance peak motion due to a magnetic field in weakly coupled chaotic quantum dots 



O 
O 

o 



oo 



c3 



-a 

o 
o 



> 

wn 
o 
o 
o 

-I— > 



i 

S3 
O 

o 



X 



IlyaL .Kurland 1 , Richard Berkovits 1,2 ' 3 and Boris L. Altshuler 1,2 

1 Physics Department, Princeton University, Princeton, NJ 08544 
2 NEC Research Institute, 4 Independence Way, Princeton, NJ 08540 
3 Minerva Center and Department of Physics, Bar-Ilan University, Ramat-Gan 52900, Israel 

(February 1, 2008, draft 1.0) 

We study the influence of moderate exchange interactions of electrons on the behavior of the peaks in the 
conductance of single electron transistors. We numerically reproduce recently observed features of both the 
peak positions and the peak heights magnetic field dependence. These features unambiguously identify the total 
spin S of each ground state. We evaluate the probability of each S (combination of spontaneous and induced 
magnetization) as a function of the exchange strength, J, and external magnetic field , B. The expressions 
involve only J and g-factor as adjustable parameters. Moreover, in a surprisingly broad parameter range these 
probabilities are determined by certain linear combinations of J and B. 

PACS numbers: 73.23.Hk, 73.23. -b,71.24.+q 



In the Pauli picture, electrons populate the orbital states of 
a chaotic quantum dot or metallic grain in a sequence of spin 
up - spin down electrons. This means that the total spin of 
the system can be either zero (when the number of electrons 
N is even) or one half (at odd 7Y). It is well known, [0 that 
measuring the electron transport through a quantum dot one 
effectively studies its ground state (GS) properties. One can 
measure, e.g., the ohmic conductance by weakly connecting 
the dot with source and drain electrodes. This conductance 
can be studied as a function of the voltage, V g , applied be- 
tween the dot and a third electrode - gate. The energy of the 
system at a given number of electrons, En, is determined by 
the gate voltage, V g . Each time when En{V 9 ) coincides with 
En+i{V 9 ) the conductance increases dramatically. Such a 
coincidence determines the position, Vg of N-th peak in 
the conductance, which is proportional to /xjy, determined as 
/j at = £jv(0) — £jv_i(0). Accordingly, the distance between 
the peaks is oc /in — /A/v-i- The height of the conductance 
peak is determined by the GS wave functions. Thus, studies of 
the conductance as function of V g allow to deduce properties 
of the manyparticle GS. 

Recently, the behavior of the conductance peaks in metallic 
dots [pp, carbon nanotubes [j^], small chaotic GaAs [||] and Si 
[||] quantum dots was monitored as a function of B. In the 
following discussion we neglect the magnetic field coupling 
to the orbital degrees of freedom (For a 2D quantum dot one 
is allowed to do it when the magnetic field is inplane. For an 
ultrasmall dot this is a good approximation for a field in any 
direction, provided that it is small enough). In other words, 
we expect the field to manifest itself only through the Zeeman 
splitting - the field shifts the GS energy by gusSB, where \ib 
is the Bohr magneton and S denotes the GS total spin. For sys- 
tems studied in Refs. the Zeeman splitting for 5=1/2 
becomes comparable to the mean single electron level spacing 



tesla for the GaAs [Q] dot and B ~ 10 tesla for the Si dot 
[§]. Thus, up to a few Tesla gy~bB <C Si, i.e., in the Pauli 
picture the magnetic field only seldomly exceeds Sj result- 
ing in changes of the GS spin. As a result one should ex- 
pect the following conductance peak behavior as function of 
the magnetic field (see Fig. |l]a): (i) The peak positions (/ij\r) 
vary with B as consecutive pairs creating an altering pattern 
of downward and upward moving peaks. The trajectories of 
the peaks are straight lines with the same slope of g/i&/2. (For 
the peak spacings (A at) the Pauli picture predicts a similar 
pattern with slopes of g/i&). Spin-orbit interactions may lead 
to fluctuations in the g factor resulting in fluctuations in the 
slope magnitudes which will not change the pattern of 
downward and upward moving peaks, (ii) The peak heights 
are the same for each pair and change between neighboring 
pairs, (iii) Once gfibB exceeds a particular spacing Sj, the 
peaks should cross. The peak height is expected to switch at 
each crossing. 



Si =< Si >; 



H+li 



(1) 



(here Si denotes the orbital energy of the one-electron orbital 
state i, and < . . . > stands for the averaging.) at B ~ 1 



Experimental results confirm this picture for metallic dots 
and carbon nanotubes [2^, but contradict it for the semicon- 
ducting chaotic dots [Q, 5 1. In Fig. |l]b we show the behavior of 
the conductance peaks as function of the magnetic field that 
follows from a model discussed later. Fig. [j]b reproduces 
features of the semiconducting chaotic dot experiments (see 
e.g., Figs. 2 and 3 in Ref. [Q) much better than Fig. |l]a. 
Peaks move as function of the magnetic field in bunches of 
two [QQj (peaks [3] and [4] in Fig. (TJb) in the same direction 
and with the same slope. For these peaks the height does not 
come in pairs. Not every change in the direction of the motion 
of a peak looks like a crossing. For example, some of these 
changes are more abrupt - a particular example is peak [3]. 
Thus, a description of the system beyond the Pauli picture is 
needed. 
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FIG. 1.: Typical conductance dependence on the gate voltage 
and magnetic field, B, with and without exchange interaction: (a) 
J = 0, (b) J = 0.15i. fi' N is determined as f/ N = /xjv — e 2 /C 
and taken in arbitrary units. The vertical curves indicate the con- 
ductance while the circles indicate the peak maxima positions. The 
behavior of specific peaks (indicated by [1] — [4]) in (b) and the dif- 
ferent regions of GS spin (depicted in the plot in units of 1/2, i.e., 1 
corresponds to S — 1/2, and separated for clarity by vertical lines) 
is discussed in the text. 

This letter aims to explain the puzzling features of the mag- 
netic field dependence of the conductance peaks positions and 
heights by spontaneous spin polarization of the electrons in 
the dot. The possibility of spontaneous magnetization has 
been in the center of much recent research It has been 

previously suggested as an explanation [^] to the puzzling ab- 
sence of a bimodal distribution in the conductance peak spac- 
ings Jj~5|]. Such a magnetization was also predicted to lead 
to kinks in the parametric motion of the peaks, (e.g., due to 
orbital effects of a perpendicular magnetic field) [O], Au- 
thors of a recent paper based on a random interaction model 
note some deviations of the peak motion from the Pauli 



picture due to spontaneous magnetization. 

The magnetization results from the exchange interaction 
between the electrons. In the presence of interaction the Pauli 
picture no longer holds. Indeed, it might be advantageous to 
lose in kinetic energy by placing electrons at higher orbitals 
in order to gain in the exchange energy. 

A symmetric dot is spontaneously polarized every time 
there is a partially filled degenerate orbital - this is the famil- 
iar Hund's rule in atomic physics. Accordingly, magnetization 
will occur in a chaotic dot once the orbitals are close enough 
to each other. Here we demonstrate that the exchange interac- 
tion explains the deviations from the Pauli picture listed above 
and also make definite predictions on the statistics of the peak 
motions which can be checked experimentally. 

We use the following Hamiltonian to describe the electrons 
in a chaotic dot [|l4]] : 



e 2 N 2 
2C 



JS{S + 1) + g^ b SB. 



(2) 



The first term in Eq. ^ is the one-partical Hamiltonian, where 
Si and rii are the energy of the i-th orbital and its occupation 
number. We assume that the energies are characterized by 
the Random Matrix spectral statistics. Here we consider both 
orthogonal (GOE) and unitary (GUE) Dyson ensembles. The 
second term represents the charging energy. We ignore fluc- 
tuations in the capacitance C JT^ ] which are not important for 
the magnetic behavior. It was demonstrated in Ref. [|l4| for 
the most general form of the electron - electron interaction 
that the sample to sample and level to level fluctuations of 
the exchange interaction are negligible for systems with high 
Thouless conductance. This is the reason why the exchange 
energy can be presented in a simple form of the third term 
of E q. (Eh with some constant J that does not fluctuate. (See 
Ref. Q| for the details.). Note that Refs. [juj], [jo]] introduce 
different exchange terms. 

We calculated the conductance using the energies En and 
wave functions obtained numerically for GOE and GUE ran- 
dom matrix realizations. An example of the peak positions 
and heights evolution as functions of the magnetic field for 
a particular GOE realization is presented in Fig [j] for (a) 
J = and (b) J = 0.1<5i. For both cases we present 
A*jv = MJV - e 2 /C. 

It is not surprising that in Fig [j]a one can see all the previ- 
ously described features of the Pauli behavior. On the other 
hand, once even a weak exchange interaction is included the 
behavior changes qualitatively. For example, even in the 
vicinity of B = 0, around a 5 — 1 state (peaks [3] and [4] 
in Fig. [j]b) the usual Pauli pattern is disrupted. Instead of an 
up and down moving pairs of identical peak height we see two 
down moving and then two up moving peaks with alternating 
height. The reason for this behavior is that first two consecu- 
tive orbitals are first filled with down spin electrons and only 
later they acquire up electrons. Generally the enhancement of 
the spin of the dot by 5 will be accompanied by a 25 bunch 
of alternating height peaks moving with the same slope. If 
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the peak spacing A at is plotted, two sets of 2S — 1 fiat curves 
sandwiching a sloped one will appear. Changes to the crossing 
patterns should be also attributed to the spontaneous magneti- 
zation. For example, such a case occurs when the spin of the 
GS between peaks [3] and [4] (indicated in Fig. [j]b) switches 
from S = 1/2 to S = 3/2. Peaks [3] and [4] abruptly change 
their slopes and amplitudes (actually the amplitude vanishes 
due to spin blockade JTo]]). The amplitude of peak [4] is equal 
to the amplitude of peaks [1] and [2] since the orbital, which 
occupation was documented by appearance of the [1] and [2] 
peaks, becomes once again available for tunneling through the 
combined effect of spontaneous and field-induced magnetiza- 
tion. Another change of slope and amplitude is seen once the 
GS between peaks [2] and [3] switches from S = to S = 1. 
Thus, the exchange interaction can explain some of the quali- 
tative behaviors seen in the experiment. 




J/5, 

FIG. 2.: Probabilities Pj(S) of S = 1 and 5" = 3/2 at B = 
as functions of the exchange strength J for GOE (/3 = 1) and GUE 
(/3 = 2). The symbols represent numerical results, the solid curves 
correspond to the first order terms in Eq. ^ while the dashed lines 
involve the higher order corrections. 



Here [-RJ denotes the integer part of a real number R and, as 
usual, p = 1(2) corresponds to GOE (GUE). In Eq. (§) we 
took into account the main asymptotics of Pj(S) at J — > 
and the first order correction in J 2 /Sf. The coefficients Cg 
and Kg depend on both (3 and S. Their values for S = 1,3/2 
are presented in table I. 

According to Eq.(|3|) the probability of a GS with a certain 
value of S depends on j3 and goes down rapidly as S increases. 
In Fig. H we present the probabilities of S = 1 and 5 = 3/2 
realizations as functions of the exchange strength, J, for the 
GOE (J3 = 1) and GUE (J3 = 2). The probabilities are ob- 
tained from a direct numerical solution of Eq. (^) (symbols) 
as well as from the analytical estimation, Eq. (|3|). One can see 
that for J < 0.3 both results coincide. 

As already mentioned, one can unambiguouslydetermine 
the spin of each GS from the trajectories of the conductance 
peaks. This fact creates an opportunity to extract J from the 
experimentally determined probabilities of different values of 
S. Note, that it is only one dimensionless parameter J/ 6\ that 
controls the probabilities Pj (S) of all spins at both possible 
values of (5. 

Magnetic field dependence of the probabilities of a given 
spin realization, Pjb(S), can be calculated in a similar way. 
It turns out that as long as the following linear combination of 
J and B, 



remains small, the probabilities behave as 



This model can also give quantitative predictions regard- 
ing the peak position movements. As illustrated above, the 
appearance of higher spin clearly manifests itself in these 
trajectories. Using Eq.(pj) one can predict the frequency of 
spontaneous magnetization appearances. For weak exchange, 
J <C Si, the probability that the dot GS has a spin S is de- 
termined by the probability of finding 25 orbitals so close to 
each other that the gain in exchange energy (due to the po- 
larization) overwhelms the loss in the kinetic energy (due to 
single rather than double occupations of the orbitals). 

Using methods described in details in [|l7|| one can show 
that a sequence of 2S levels with a total energy smaller than 
JS(S + 1) which is necessary for the spin to reach the value 
S at J <C Si takes place with the probability: 

/ 7 \ (/3S+1)(2S-1) , rf\ 



P JB (S) = G f jrO»s+i)(2S-i) (1 _ K X ^), (5) 



In Fig. |a the probabilities Pj, B (5) for S = 1, 3/2 calculated 
numerically at J/ Si = 0.1 are compared with the asymptotics 
of Eq. (^|). The agreement turns out to be good as long as 
gfj, B BS < 0.4$ x . 

The form of Eq. (||) and the higher order terms suggests that 
exchange interaction and Zeeman splitting enter the probabili- 
ties Pj t s (S) only through the combination X. As can be seen 
from Fig. ||b, a general scaling 



Pj,b(S) = F(S,X), (6) 
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FIG. 3.: (a) Probabilities of S = 1, 3/2 for J = 0.1<5i as functions 
of the magnetic field B. The symbols represent numerical results, 
while the curves (solid - first order, dashed - second order) represent 
Eq. (Q). (b) The same probabilities of measuring a spin value S for 
different values of B and J scaled as functions of the parameter X, 
(Eq. (Q)). Hollow symbols represent the GOE and filled symbols 
stand for the GUE. Large (small) symbols correspond to S = 1 
(S = 3/2). One can clearly see that the scaling form holds much 
beyond the first couple of terms estimated here. Inset: probability 
of S = 1, in a /3 = 1 Hubbard model. The Hubbard model was 
solve using exact diagonalization (similar to Ref. [^]) for different 
lattices and fillings (4 x 3, 4 x 4, 4 x 5 with 4, 6, 8 electrons ), 
and different values of on-site interactions, U — — 2t, and mag- 
netic fields. The scaling parameter X was calculated using J(U) 
deduced from the spin dependence of the averaged energy. The x 
symbols correspond to the Hubbard numerical results, the y are the 
same as in (b). A full description of the Hubbard model results will 
be given elsewhere. 

where F(S,X) is some function of the scaling parameter 
X and the spin S, holds for values of X which much exceeds 
the range of validity of Eq. (§). Deviations from the scal- 
ing law Eq. (^]) at higher X are probably due to the fact that 



the probability to observe even higher spins (S > 3/2) be- 
comes already substantial. At J close to 6% one should use the 
approach of Ref. Jl4] ] to describe Pj.b{S) analytically. It is 
important to stress that this scaling will hold for any system 
of interacting electrons as long as the Thouless energy is large 
and there is no spin orbit scattering. This is demonstrated in 
the inset of Fig [|b where results for a Hubbard model are pre- 
sented. 

In conclusion, we have shown that the presence of even 
a relatively weak exchange interaction explains many of the 
deviations from the Pauli picture seen in recent experiments 
[^,|]| . We calculate how the frequency at which different val- 
ues of spin appear depends on weak exchange interaction and 
a moderate magnetic field. The strength of the exchange inter- 
action, J, and the effective ^-factor are the only adjustable pa- 
rameters that determine the probabilities for a GS of the dot to 
have any given spin at any given magnetic field for both GOE 
and GUE cases. In particular we predict that these probabili- 
ties will follow a one parameter scaling law over a wide range 
of magnetic fields and exchange interaction strengths. To 
test these predictions experimentally one needs only enough 
statistics for the behavior of the conductance peaks. 
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S = 1 


S = 3/2 


S= 1 


5 = 3/2 


c 

K 


tt73 

7T 2 /5 


97r 4 /50 
18tt 2 /49 


8tt79 
8tt 2 /25 


81ir e /400 
792tt 2 /1225 



TABLE I.: The factors Cg and Kg appearing in Eq. (^). 



